Charmonium potentials: Matching perturbative and lattice QCD 
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Central and spin-spin potentials for charmonium, constructed from Nambu-Bethe-Salpeter am- 
plitudes in lattice simulations of full QCD, are matched with results from perturbative QCD at an 
appropriate distance scale. This matching is made possible by defining the perturbative potentials 
through Fourier transforms with a low-momentum cutoff. The central (spin-independent) poten- 
tial is compared with potentials derived from an expansion in powers of the inverse quark mass. 
A well-controlled continuation of the charmonium spin-spin potential from lattice QCD to short 
distances is performed. The mass splittings of the charmonium singlet and triplet states below the 
open charm threshold, obtained from the matched spin-spin potential, are in good agreement with 
the experimental values. 



I. INTRODUCTION 

Within the framework of potential non-relativistic 
QCD (pNRQCD) the heavy-quark-antiquark potential is 
a well-defined quantity, and usually presented as an ex- 
pansion in powers of the inverse heavy-quark mass 
The leading-order static potential (i.e., the potential be- 
tween infinitely heavy quarks) is accurately known from 
lattice QCD studies. The potentials at order 1/m and 
1/to 2 have been computed in quenched lattice QCD us- 
ing the Wilson- loop formalism In this approach, 
a spin-spin potential is found at order 1 /to 2 with a sign 
such that the mass ordering of hyperfine multiplets is 
opposite to the ordering observed empirically. An earlier 
extraction of the spin-spin potential from lattice QCD [f| 
found a sign in agreement with the empirical ordering. 

A different lattice QCD approach to extract the heavy- 
quark potentials has recently been proposed in Ref. [6J. 
The spin-spin and central potentials at finite quark mass 
are obtained from Nambu-Bethe-Salpeter (NBS) ampli- 
tudes through an effective Schrodinger equation. A re- 
pulsive spin-spin potential which shifts spin triplet states 
upward and spin singlet states downward is found as re- 
quired for reproducing the experimental charmonium and 
bottomonium spectra. Although some model dependence 
is involved in this approach, the resulting shape of the po- 
tential is quite different from a (^-function potential that 
one commonly obtains from one-gluon exchange with an 
effective coupling strength a s treated as an adjustable 
parameter. 

In this Letter we show that one-gluon exchange is, 
in fact, sufficient to derive a perturbative spin-spin po- 
tential, consistent with the recent findings in lattice 
QCD, if the running of the QCD coupling a s (q) is prop- 
erly included. A matching of perturbative and non- 
perturbative regions of the spin-spin potential is thus 
possible. The central (spin-independent) charmonium 
potential, obtained in full lattice QCD from NBS ampli- 
tudes, is matched to a corresponding perturbative poten- 
tial and compared with results obtained previously in the 
1/m expansion JJ. The perturbative potentials are con- 
structed with a restricted Fourier transformation. This 
method replaces the frequently used renormalon subtrac- 



tion scheme 0] . The S-wave charmonium spectrum is de- 
rived in Section [IIII from the combination of the matched 
central and spin-spin potentials and compared with ex- 
perimental values. The single free parameter in our ap- 
proach is an overall additive constant which enters in the 
extraction of the charm-quark mass. 



II. CONSTRUCTION AND COMPARISON OF 
CHARMONIUM POTENTIALS 

The quark- antiquark potential is commonly written in 
the form 



V qq (r) = V c (r) + S q -SgV s (r) + 



(1) 



with a central potential Vc{r) and a spin-spin poten- 
tial Vs(r). First, consider the perturbative central (spin- 
independent) charmonium potential. As in Ref. [7|, we 
define it in r-space via a restricted Fourier transformation 
with a low-momentum cutoff /ic'- 
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with q — \q\. It includes the static potential V^(q) at 
two-loop order in momentum space 
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The coefficients a\ and ai in Eq. ([3]) are known analyti- 
cally [lfj| - [T2l | and have the values 
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FIG. 1. Charmonium potential from a combination of pertur- 
bative QCD and lattice QCD [16( matched at r m = 0.14 fm 
(solid curve). The dot-dashed line (with error band) shows 
the static potential from Ref. 0] for comparison, while the 
dashed line (with error band) shows the static-plus- 1/m po- 
tential from Ref. 17]. The energy scale is chosen relative to 
the potential at r — 0.5 fm for convenience. 



for three light quark flavors. The low-momentum region 
excluded in Eq. is not accessible in perturbation the- 
ory. It is substituted by an additive overall constant to 
the potential (see Q] for details). 

The full four-loop renormalization-group (RG) run- 
ning [l3| with three light flavors is implemented for 
the strong coupling a s (q). We use the input value 
a s (1.25 GeV) = 0.406 ± 0.010, derived from a s (m z = 
91.1876 GeV) = 0.1184±0.0007 [3, taking into account 
flavor thresholds fljj . 

The perturbative potential is matched at a suitable dis- 
tance r m to the spin-independent potential calculated in 
full lattice QCD using NBS amplitudes 0. This lattice 
potential includes charm-quark mass effects to all orders 
and can be parametrized as 



V' at (r) 



A 



a r + const . 



(7) 



with A = 0.813 ± 0.022 and = (0.394 ± 0.007) GeV. 
Optimal matching is achieved with a low-momentum cut- 
off nc = (0-54 ± 0.02) GeV. In the region around the 
chosen matching position r m — 0.14 fm, the perturbative 
potential and the lattice potential are both expected to 
be reliable. The potential changes only marginally under 
limited variations of the matching point. 

In Fig. Q] we compare this matched central potential 
(solid curve) with the potential obtained in the Wilson- 
loop approach using a 1 /m expansion. One observes that 
the "full" potential differs significantly from the static po- 
tential constructed in Ref. 7] (dot-dashed line with error 
band), but it is consistent within errors when the 1/m 
potential from Ref. [7| is added to the static potential 
(dashed line with error band). The charm-quark mass m 
relevant for the 1/m potential has been varied in Fig. Q] 
in the range (1.5 ± 0.2) GeV. 



Let us now focus on the spin-spin potential. Tensor 
and spin-orbit terms are not discussed in this Letter. 
These have so far not been studied in the new lattice 
QCD approach employing the NBS amplitudes. The au- 
thors of Ref. [IH fit the spin-spin potential from full lat- 
tice QCD with three different functional forms. We use 
the exponential form 
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with a = (0.825 ± 0.019) GeV and £ = (1.982 ± 
0.024) GeV, since it provides the best fit to the lattice 
data. 

In the following we construct a perturbative spin-spin 
potential that can be used to continue the lattice po- 
tential to short distances. Recall first the well-known 
(schematic) spin-spin potential derived from one-gluon 
exchange assuming a constant coupling a s : 
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with the quark/ antiquark mass m q . This <5-function term 
is in agreement with the leading-order spin-spin potential 
obtained in pNRQCD [J. 

In analogy to the case of the spin-independent poten- 
tial, it is possible to include the RG running of a s (q) in 
the construction of the spin-spin potential. We define the 
perturbative part of Vs(r) as: 
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with a low- momentum cutoff us- For q > fi$ the pertur- 
bative RG evolution of a s (q) is supposed to be reliable. 
The non-perturbative infrared behavior of the quark and 
gluon couplings prohibits a controlled low-momentum ex- 
tension for q < us. It is nevertheless useful to examine 
such an extension for r <C i n the form 
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with a parameter a s reflecting the average interaction 
strength in the infrared region. Note that Eq. (|11[) gives 
a positive constant proportional to to be added as a 
correction to Vg ert (r, /is), with a s not known, but ex- 
pected to be of 0(1). 

In order to facilitate the numerical evaluation of the 
Fourier integral for the perturbative part, it is useful to 
rewrite Vj >crt (r, /is) in the form 



V s pert (r, 
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To be consistent with the lattice QCD analysis, the 
charm-quark mass m in the denominator of Eq. (|12l) is 
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FIG. 2. Spin-spin potential in r-space derived from one-gluon 
exchange as denned in Eq. <|12[) . The potential is shown for 
different values of the low-momentum cutoff /is. 
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FIG. 3. Charmonium spin-spin potential (with error band) 
from a combination of perturbative QCD (see Eq. ([12}) and 
lattice QCD [16|], matched at r m = 0.14 fm. Dashed line: con- 
tinuation of the exponential lattice fit ([8} to short distances. 



identified with (1.74 ±0.03) GeV, the kinetic quark mass 
that has been determined in Ref. (l6| . Discarding the 
additive constant Vg in the first step, the resulting short 
distance potential and its dependence on the cutoff scale 
/is is shown in Fig. [5J The shape is evidently very dif- 
ferent from a Gaussian or a 5-function. Such forms are 
frequently used for the spin-spin potential in phenomeno- 
logical models. 

The two components of the spin-spin potential, arising 
from perturbative QCD for r < r rn and from lattice QCD 
for r > r m , can be matched at r m = 0.14 fm. When ig- 
noring the additive constant Vg this matching is achieved 
for /i S = (0.43 ± 0.02) GeV. With inclusion of the (vari- 
able) additive constant V s r the infrared cutoff /is can be 
varied over a wide range almost without any effect on the 
spin-spin potential. Figure [3] shows the matched poten- 
tial with an infrared cutoff of /is = (0.75 ± 0.25) GeV. 



III. SPECTROSCOPY AND CHARM-QUARK 
MASS 

Here, we discuss the charmonium spectrum below the 
DD threshold as derived from the matched potentials 
constructed in the previous section. We focus on the 
IS and 2S states which are not influenced by the tensor 
and spin-orbit interactions. The Schrodinger equation 
for these states, 



1 /m-correction term: 
mpsivc) = ™ P oic 

+ 1 /^U)( g) + ^M1 (i4) 

2 J (2tt) 3 [ W m/2 J 1 ' 

q<Hc 

Replacing the operator S q -Sq by its eigenvalues, —3/4 
for the spin singlet and 1/4 for the spin triplet, the 
Schrodinger equation (|13j) is solved numerically. The 
spin-spin potential Vs(r) can be included in two different 
ways. In the first case it is treated in first-order pertur- 
bation theory, in the second case it is fully included in 
the Schrodinger equation. Due to the singular behavior 
of Vs(r) ~ r ~ 2 - 8 f or r _ > o (according to our construc- 
tion) the wave functions diverge (mildly) in the latter 
case for very small values of r. In this case we solve the 
radial Schrodinger equation numerically for r > 0.003 fm 
and convince ourselves that the physical results are not 
affected by contributions coming from shorter distances. 
The results of the two alternative treatments of the spin- 
spin potential agree within 12 MeV (see Table fl}. The cal- 
culated mass splittings between singlet and triplet states 
are in good agreement with experimental results for both 
IS and 2S charmonia. 

The single free parameter mps{nc) is chosen such that 
the spin- weighted average of the IS states agrees with its 



V 2 
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involves a single free parameter mps(/ic), the (nc~ 
dependent) charm-quark mass in the potential sub- 
tracted (PS) scheme [T7| . In order to make the PS scheme 
applicable for our purposes, we extend it by including the 









Case 1 


Case 2 


Experiment [18] 


IS 


mass splitting 


[MeV] 


117 ±6 


105 ±6 


116.6 ± 1.2 


2S 


mass splitting 


[MeV] 


56 ±3 


46 ±3 


49 ±4 



TABLE I. Predicted mass splittings of charmonium IS and 
2S multiplets in comparison with experimental data. Case 1: 
spin-spin potential treated in first-order perturbation theory. 
Case 2: spin-spin potential fully included in the Schrodinger 
equation. 
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FIG. 4. Predicted masses of charmonium IS and 2S states 
in comparison with experimental data [18|]. The results in 
the first column are based on the spin-independent potential 
only. The effects of the spin-spin potential, treated in first- 
order perturbation theory, are added in the second column. 



experimental value {T8|. For the excited 2S states our 
approach predicts masses for ?7 C (2S) and ^>(2S) that are 
slightly too large (see Fig. |4|). However, these states are 
close to the DD threshold. Going beyond this threshold 
requires a complex (energy-dependent) cc potential or an 
explicit treatment of coupled channels. While the imag- 
inary part starts at the opening of the DD channel, the 
corresponding dispersive real part induces an attractive 
shift of the 2S states. In second-order perturbation the- 
ory this shift is proportional to the squared cc — > DD 
transition matrix element. 

The predicted size of the hyperfine splittings is quite 
sensitive to the value of the mass to in the denominator 
of Eq. (jT2"j) . For example, choosing to = 1.5 GeV instead 
of the kinetic quark mass (1.74 ± 0.03) GeV [16| would 
give rise to a IS mass splitting that is about 20% too 
large. This is in contrast to variations of the matching 
position r m and the infrared cutoff [is, which affect the 
hyperfine splittings only marginally. 

The spin-spin potential, as constructed in the previ- 
ous section, produces a non- vanishing but small splitting 
between the IP singlet and triplet states, namely h c (lP) 
and Xcj(lP), unlike the ^-function spin-spin potential. 
In first-order perturbation theory the effect amounts to 
a mass difference of (8.2 ± 0.5) MeV. The full inclu- 
sion of the spin-spin potential Vs(r) in the Schrodinger 
equation gives rise to a slightly larger mass splitting of 
(8.3 ± 0.5) MeV. 

The value of the mass parameter mps(nc), determined 
in our approach by fitting to empirical charmonium spec- 
tra, can be translated into alternative schemes for quark 
masses. The PS mass mps(pc) is first converted to the 
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This procedure is described in detail 
in Ref. 7]. Applying the same method here, we find for 
the charm-quark mass in the MS scheme 



rn c 



(1.21 ±0.04) GeV, 



(15) 



in good agreement with other determinations 

UM- The 

error reflects combined uncertainties in the lattice poten- 
tials, in the input value of the strong coupling a s {q), and 
from the matching to the empirical states. 

We close with a few remarks concerning bottomo- 
nium: until now, the bottomonium spin-spin potential 
has not been studied within the new lattice QCD ap- 
proach based on NBS amplitudes. An extrapolation of 
the spin-spin potential from charmonium to bottomo- 
nium can be done by simply assuming a 1/m 2 depen- 
dence of the lattice potential and allowing for varia- 
tions of the mass parameter to. In the perturbative part 
of the potential we account furthermore for a modified 
running of ot s {q) due to four massless flavors and use 
a s (4.2 GeV) = 0.226 ± 0.003 as an input. The empirical 
mass splitting of (69 ± 3) MeV pj] between ryb(lS) and 
T(1S) can be reproduced either for a kinetic bottom- 
quark mass m = 4.7 GeV (with the spin-spin poten- 
tial treated in first-order perturbation theory), or with 
to = 4.3 GeV (if the spin-spin potential is fully included 
in the Schrodinger equation). It will be interesting to 
have available the corresponding lattice QCD results for 
bottomonium. 



IV. SUMMARY 

Central and spin-spin potentials for charmonium have 
been derived by combining perturbative QCD at small 
distances (r < 0.14 fm) with results from lattice QCD 
for larger distances up to r ~ 1 fm. By defining the 
perturbative potentials via a restricted Fourier transfor- 
mation this matching has been made possible. We have 
found that the central quark-antiquark potential, con- 
structed from NBS amplitudes in full QCD lattice simu- 
lations [f| [l6j , agrees within errors with the static-plus- 
\Jra potential derived in the Wilson-loop formalism 0- 
4]. The matched spin-spin potential produces hyperfine 
splittings for the S-wave charmonium states that are in 
good agreement with experiment. The MS mass of the 
charm quark also agrees well with other determinations 
of this QCD parameter. 
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